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The present work reports study on the interacting Ricci dark energy in a modified gravity 
theory named f(R,T) gravity. The specific model f(R,T) — jiR + vT (proposed by R. 



Myrzakulov, larXiv:120 5.5266v2 



29() is considered here. For this model we have observed a 



quintom-like behavior of the equation of state (EoS) parameter and a transition from matter 
dominated to dark energy density has been observed through fraction density evolution. The 



a 

60 of the universe. 



statcfinder parameters reveal that the model interpolates between dust and ACDM phases 



^ ■ I. INTRODUCTION 

nri nn 

The origin of dark energy (for review see [1H2J) responsible for the cosmic acceleration [4|, [5J is 
one of the most serious problems in modern cosmology. The first step toward understanding the 
nature of dark energy is to clarify whether it is a simple cosmological constant or it originates from 
other sources that dynamically change in time [6]. In an extensive review, Nojiri and Odintsov 
[3] thoroughly discussed the reasons why modified gravity approach is extremely attractive in the 
applications for late accelerating universe and dark energy. Other remarkable reviews on modified 
6|, i Various modified gravity theories have been proposed so far. These include, f(R) 

theories. 



gravity are 
3,0, f(T) 
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14], f(G) 0,0, Horava-Lifshitz 0,0 and Gauss-Bonnet 0, 
The present work concentrates on f(R,T) gravity, with T being the trace of stress-energy tensor, 
manifesting a coupling between matter and geometry. Before going into the details of f(R,T) 
gravity, let us first briefly survey the f(R) gravity. The recent motivation for studying f(R) 
gravity has come from the necessity to explain the apparent late-time accelerating expansion of 



the Universe [8j. Some extensive reviews of f(R) gravity are 2ll424||. Thermodynamic aspects of 
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25 



261 ] . A generalization of /(-R) modified 



/(i?) gravity have been investigated in the works of 
theories of gravity including in the theory an explicit coupling of an arbitrary function of the Ricci 
scalar R with the matter Lagrangian density C m leads to the motion of the massive particles is 
non-geodesic, and an extra force, orthogonal to the four- velocity, arises 27], [28[. Harko et al 27] 
proposed an extension of standard general relativity, where the gravitational Lagrangian is given 
by an arbitrary function of the Ricci scalar R and of the trace of the stress-energy tensor T and 
dubbed this as f(R, T) gravity. The f(R, T) gravity model depends on a source term, representing 
the variation of the matter stress-energy tensor with respect to the metric. A general expression 
for this source term is obtained as a function of the matter Lagrangian £ 
paper, Myrzakulov 



271 ] . In a recent 



291 ] derived exact solutions for a specific model f(R, T) = fiR + vT and showed 
that for some values of the parameters the expansion of our universe can be accelerated without 
introducing any dark component. The present work aims to reconstruct the Ricci dark energy 



(RDE) 



30 



32 



34 



381 ] under f(R, T) gravity. Rest of the work is organized as follows: In section 



II we have briefly reviewed RDE. In section III we have presented an overview of f(R,T) gravity. 
In section IV we have reconstructed interacting RDE in f(R,T) gravity. We have concluded in 
section V. 



II. A BRIEF OVERVIEW OF RICCI DARK ENERGY 

Gao el.al 3(| proposed a holographic dark energy model in which the future event horizon is 
replaced by the inverse of the Ricci scalar curvature, and dubbed this model the "Ricci dark energy 
model" (RDE). This model (i) avoids the causality problem (ii) is phenomenologically viable, and 
(iii) can solve the coincidence problem of dark energy 32| . The Ricci curvature of FRW universe 
is given by (32 ] 

k 



R 



6 { H + 2H + 



(1) 



where, k is the curvature of the universe and a is the scale factor. The energy density of RDE is 
given by 



35] 



Px 



In flat FRW universe, k = and hence we have 



3c 2 ( H + 2H 2 + 4 



(2) 



px = 3c 2 [ H + 2H 2 



(3) 
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In the present work we are considering RDE interacting with pressureless dark matter with energy 
density p m . Various forms of "interacting" dark energy models have been constructed in order to 
fulfil the observational requirements. Plethora of literatures are available where the interacting 
dark energies have been discussed. Several examples of interacting dark energy are presented in 



3614411] . In a subsequent section we shall consider the interacting RDE in f(R,T) gravity. The 



metric of a spatially flat homogeneous and isotropic universe in FRW model is given by 

ds 2 = dt 2 - a 2 (t) [dr 2 + r 2 {d6 2 + sin 2 fl# 2 )] (4) 
where ait) is the scale factor. The Einstein field equations are given by 

H 2 = \p (5) 

and 

H = -\{p + p) (6) 

where p and p are energy density and isotropic pressure respectively (choosing 8irG = c = 1). The 
conservation equation is given by 

p + 3H{p + p) = (7) 
As we are considering interaction between RDE and dark matter, 

P = Px + Pm, P = Px (8) 

It should be stated that we are considering pressureless dark matter, p m = 0. Since the components 
do not satisfy the conservation equation separately in the case of interaction, we need to reconstruct 
the conservation equation by introducing an interaction term Q. The interaction term could be in 



any of the forms 42]: Q oc Hpx, Q oc Hp m and Q oc H(px + Pm)- In the present paper we take the 
interaction term in the second of the three forms mentioned above. Accordingly the conservation 
equation is reconstructed as 

px + 3H{p x + px) = 3H5p m (9) 



p m + 3Hp m = -3H5p rn 



(10) 



III. THE f{R,T) MODEL 



One of interesting models of f(R, T) gravity is the so-called M37-model. Its action is 29j 



S = J f(R,T)^d 4 x + J CmV^gd^x 



(11) 



where C m is the matter Lagrangian and f(R,T) is an arbitrary function of R and T, where R is 
the scalar curvature and T is the torsion scalar. Here, 



R = u + 6(H + 2H 2 



(12) 



T = v-6H 2 



Ref [29(] considered the following model of f(R,T) 



f(R, T) =nR + vT 



(13) 



(14) 



where \jl and v are real constants and u and v are taken as u = aa n and v = (3a m with m, n, a 
and f3 as real constants. The equations system of this f(R,T) is 



liDi + vE x + K(uT + /iR) = -2a 3 p 



(15) 



pAx + vBx + M(vT + (iR) = 6a 2 p 



(16) 



where, 



p + 3H(p + p) = 



D\ = a 3 ( 6 — h au a 



. a 
a 



(17) 



(18) 



E 1 =a A \ -12^r + av c 



K = -a 6 



A\ = 12a 2 + 6aa + 3a 2 aua + a 3 ^ — a 3 ii a 



i?i = — 24d 2 — 12aa + 3a 2 ava + a 3 Va — a 3 v a 



(19) 
(20) 
(21) 
(22) 
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M = -3a 2 



(23) 



R = u + 6(H + 2H 2 



(24) 



T = v — 6H 



Subsequently, the modified field equations are obtained as [29 ] 

3(/j + v)H 2 + ]^{paa n + vfia m ) = p 



( M + u)(2H + 3F 2 ) + Mn + 3) q" + + 3) 



6 



6 



a m = -p 



(25) 

(26) 
(27) 



IV. INTERACTING RDE IN THE f(R,T) GRAVITY 



As stated earlier, p = px + Pm and p = px are taken in the equations (26) and (27). From 
equation (10) we get 



Pm — Pmoa 



-3(1+5) 



(28) 



Using equations (3) and (28) in the right hand side of the equation (26) we get the Hubble's 
parameter as a function of the scale factor a as 

3+m (3va 3+n 2a- 35 p m0 



2(-2c 2 + M + 1 /) 1 



apa 



+ 



+ 



c 2 (4 + m) -2{ji + v) c 2 (4 + n) -2(p + u) c 2 (-l + 35) + 2(p + i/) 

(29) 



Subsequently we get and H as functions of a as follows 



H 



^ d(-2c 2 +ii+u) 
-71 



a c 



, a m man , a n n$v -3(i+a)(i+f) Pm0 

~ l "6(c 2 (4+m)-2(^+;y)) 6(c 2 (4+n)-2(^t+z/)) c 2 (-l+35)+2(n+v) 



(30) 



2(m+^) 

12a c 2 Ci(-2c 2 +/j+z/) 2 



, ■■ , a 4+n n 2 l3u , ISa 1 - 36 (1+8) 2 p m0 

_l "c 2 (4+m)-2(/ t +5) + c 2 (4+n)-2(At+^) + c 2 (-l+35)+2(/x+*/) 



(31) 



Using (29) and (30) in equation (3) we get the energy density for RDE under interaction with 
pressureless dark matter under f(R,T) gravity as 



1 



Px 



2(-2c 2 + M +,) a m c 2 apU + m) a n c 2 f3vU + n) 2a- 3 ( 1+<5 ) C 2 (-l + 35)p m0 

OO + V) + 0/J ; — ; ^ + 



c 2 (4 + m) - 2(p + u) c 2 (4 + n) -2{p + u) c 2 (-l + 35) + 2{p + v) 

(32) 
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Using the above form of px in the conservation equation (9) we get the pressure for RDE in the 
present case as 



Px 



-a n (3 + n)ap - a m (3 + m)$v - 6(/i + v) 



a m (3+m)a/x 



3(c 2 (4+m)-2((U+y)) + Z(c 2 {i+n)-2(n+u)) c 2 (-l+38)+2(n+u) 



a"(3+n)ffy 



2(-2c / + M + i/) 

73 r- 



(33) 



Using the expressions for energy density and pressure derived above we get the equation of state 
parameters 

Px 



and 



w x 



Px 



(34) 



W total 



PX 



PX + Pm 



; (Pm = o) 



(35) 



* -0.8 




FIG. 1: Behavior of wx = ^ against redshift z 
a- 1 -!. 



The deceleration parameter q 



43J comes out to be 



aa 



H 
IP 



S.2 



where, 

6a 2 (- 2c2 +^+^Ci(-2c 2 + /i + v) 

Ci = -3 + 



a m map 



-+- 



a n nfiv 



(36) 



6a -3(l+<5) (1 + g)/Qm0 

c 2 (4 + m) - 2(/i + i/) ' c 2 (4 + n) - 2(/z + i/) c 2 (-l + 35) + 2(/i + v) 

(37) 



C2 = 2 3a" 



2(-2e' ! + M + iy) 



"Ci + 



a ma ix 



a n /3u 



+ 



c 2 {A + m) -2{p + v) c 2 (4 + n) - 2{p + v) c 2 (-l + 36") + 2(p + i/) 

(38) 
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L 



FIG. 2: Behavior of w total = 



Next, we consider the statefinder parameters {r, s} 
(30) and (31) we get the statefinder parameters as 



44j for the present case. Using equations (29), 



where, 



H H 0i 



(39) 



£•1 



6a g Ci{2e 4 -5c 2 (>+^)+2(/j+^) 2 } 



a m m(3+m)a/j a"n(3+n)/3i/ 18q- 3 ( 1 +' 5 ) (l+5)5p m0 



(c 2 (4+m)-2(p+^)) + (c^(4+n)-2(/i+i/)) c 2 (-l+35)+2(/i+;/) 



(40) 



02 = 2 



3Cia ? + 



^ { 



apa 1 



3 + m 



c 2 (4+m)-2(/i+i/) 



+ 



(Sua 



3+n 



c 2 (A+n)-2(fi+v) + c 2 (-l+3<5)+2(^+ 



2a M p m o 



(41) 



3F# + H Ci 



3if(2# + 3# 2 ) 6 



(42) 



6 



a m m(3+m.)oA t 



6(c 2 (4+m)-2(/i+^)) + 6(c 2 (4+n)~2(^+!/)) 



a n n(3+n)/3i^ 



2(-2c' i +/i+i/) 
1 " ^ 



+4a ? Ci(2c 4 - 5c 2 (^ + i/) + 2(// + i/) 2 ) + 



c 2 (-l+35)+2(p+;/) 



(43) 



. 2(-2e^+ M +i^) 

6 = 3 ( -£a 3 



a"(3+ra)/3i/ 



+ 



a m (3+m)a/j 
3(c 2 (4+ra)-2(>+i/)) 



2a- 3 ( 1 +^gp m0 



3(c 2 (4+n)-2(/i+i/)) c 2 (-l+3<5)+2(ji+i/) 



(44) 
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V. DISCUSSIONS 



In figure 1 we have presented the EoS parameter wx = §^ for RDE under f(R,T) gravity 
against redshift z = a -1 — 1. In this and the subsequent figures the solid, dashed and dotted lines 
would correspond to c 2 <, =, > 0.5 respectively. Figure 1 shows that for all values of c 2 the 
EoS parameter transits from wx > — 1 to wx < — 1 i.e. from quintessence to phantom. From 
this figure we see that at early times, roughly z > 2, the EoS approaches 0; i.e., in this model the 
dark energy behaves like dust matter during most of the epoch of matter domination. The EoS 
crosses phantom crossing wx = — 1 at z ~ —0.2 and in the distant future, the equation of state 
approaches wx = —1-2, the Universe evolves into the phantom-dominated epoch. For this model, 
the EoS crosses —1 , so it may be classified as a "quintom". Thus, the interacting RDE behaves like 



quintom in the f(R, T) gravity model proposed by [29(. In figure 2 we have plotted 



Ptotal 



Px 



PX+Pn 



where we found similar crossing of the phantom divide Wtotal = — 1 and transition from w total > — 1 
at higher redshift to Wtotal < — 1 at lower redshifts. It might be stated that we have chosen the 
model parameters as a = 1.2, (3 = 1.2, v = 0.3, \x = 0.4, n = 3, C\ = 3.02, 5 = 0.05, p m o = 0.23 
and m = 2. In all the figures, the solid, dashed and dotted lines correspond to c 2 < 0.5, = 0.5 
and > 0.5 respectively. 
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FIG. 3: Behavior of fractional densities Vtx = (red 
lines) and Sl m — (blue lines) with evolution of the 
universe. 



In figure 3 we have plotted fractional densities Ox = an d 



against redshift z. 



where, H 2 = (/i + v)H 2 + |(/xq(1 + z)~ n + uj3(l + z)~ m ). We observe that at from higher to 
lower redshifts the fractional density fix of RDE is increasing, while the fractional density of dark 
matter is decreasing. This indicates the transition from the matter dominated to dark energy 
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dominated universe. At very early stage of universe z > 2, the fractional density of dark energy 
£lx is dominated by fractional density of dark matter £Idm- After z = 2, the £lx starts showing 
an increasing pattern and £Idm starts showing a decaying pattern. This indicates the gradual 
transition from matter dominated era to the dark energy dominated era. We denote the cross-over 
point of the fractional densities by z cross , where tlx = &>dm- For c 2 < 0.5, = 0.5 and > 0.5 the 
Zcross ~ 0.6, 0.75 and 0.5 respectively. It is also observed that in the early universe the density 
contribution of dark energy can occupy roughly 20%-30%. However, at this stage the dark energy 
behaves like dust matter. So, effectively speaking, the matter density contribution is still 100%. 
Finally, from figure 3 our observation is that RDE in f(R,T) gravity is capable of achieving the 
transition from matter-dominated to dark energy-dominated universe. 
In figure 4 we have plotted the deceleration parameter q as a function of the redshift z. We 



o.o 



-0.5 



-1.0 - 



-1.5 - 




FIG. 4: Behavior of deceleration parameter q (Eq. 36). 

observe that at very early stage, roughly z > 2, q > i.e. the decelerated universe. At z « 2.5 
the deceleration parameter transits from positive to negative level. That is, the universe gradually 
transits from decelerated to accelerated stage. At later stage q = —1.5. Thus, we observe that it 
is possible to achieve the accelerated phase of the universe from decelerated phase for RDE under 
f(R,T) gravity. 

In figure 5 we have created the r — s trajectories for the three values of c 2 under consideration. 
Sahni et al. [44] demonstrated that the statefinder diagnostic could effectively discriminate different 
forms of dark energy. Cosmological constant, quintessence, Chaplygin gas, and braneworld models 
were investigated by [45J using the statefinder diagnostics and it was observed that the statefinder 
pair could differentiate between these models. An investigation on statefinder parameters for 



differentiating between dark energy and modified gravity was carried out in 



diagnostics for f(T) gravity has been investigated in [43]. In the {r, s} plane, s > corresponds 



46J. Statefinder 
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FIG. 5: The r- s trajectory (Eqs. 39 and 42). 



to a quintessence like dark energy and s < corresponds to a phantom-like dark energy, and an 
evolution from phantom to quintessence or inverse is given by a crossing of the point (r = 1, s = 0) 



in {r,s} plane 



471 ] . Also, the fixed point (r = 1, s = 0) corresponds to ACDM scenario. In figure 



5 we clearly observe a transition from quintessence to phantom as the r — s trajectory transits 
from positive to negative sides of s after crossing the (r = l,s = 0) point. Also, we find that, for 
finite r, s — > — oo that corresponds to dust phase. Thus, the interacting RDE in /(, T) gravity 
with f(R,T) = fiR + vT interpolates between dust and ACDM phases of the universe. Also, the 
statefinder diagnostics supports the quintom-like behavior of the equation of state. 



VI. CONCLUDING REMARKS 



In this work we considered interacting Ricci dark energy in f(R,T) = /iR + vT gravity. We 
have observed that the EoS parameter exhibits quintom like behavior for this model. Also, the 
transition from matter dominated to dark energy dominated universe is achievable by this model. 
The deceleration parameter have exhibited a transition from positive to negative sign, thereby 
showing the evolution of the universe from deceleration to acceleration. The statefinder diagnostics 
have been investigated and an interpolation between dust and ACDM phases of the universe has 
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been observed under this model. 
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